Nuclear Corrections in Neutrino-Nucleus Deep Inelastic Scattering and their 
Compatibility with Global Nuclear Parton-Distribution-Functions Analyses 
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We perform a global x^-analysis of nuclear parton distribution functions using data from charged 
current neutrino-nucleus {uA) deep-inelastic scattering (DIS), charged-lepton-nucleus {£^A) DIS, 
and the Drell-Yan (DY) process. We show that the nuclear corrections in vA DIS are not compatible 
with the predictions derived from £^A DIS and DY data. We quantify this result using a hypothesis- 
testing criterion based on the distribution which we apply to the total sts well as to the of 
the individual data sets. We find that it is not possible to accommodate the data from uA and £^A 
DIS by an acceptable combined fit. Our result has strong implications for the extraction of both 
nuclear and proton parton distribution functions using combined neutrino and charged-lepton data 
sets. 

PACS numbers: 12.38.-t,13.15.+g,13.60.-r,24.85.+p 



High statistics neutrino deep-inelastic scattering (DIS) 
experiments have generated significant interest in the lit- 
erature as they provide crucial information for global fits 
of parton distribution functions (PDFs). The neutrino 
DIS data provide the most stringent constraints on the 
strange quark distribution in the proton, and allow for 
flavor decomposition of the PDFs which is essential for 
precise predictions of the benchmark gauge boson pro- 
duction processes at the LHC. Moreover, the neutrino 
experiments have been used to make precision tests of 
the Standard Model (SM). A prominent example is the 
extraction of the weak mixing angle Ow in a Paschos- 
Wolfenstein type analysis A good knowledge of the 
neutrino DIS cross sections is also very important for 
long baseline experiments of the next generation which 
aim at measuring small parameters of the mixing ma- 
trix such as the mixing angle ^13 and eventually the CP 
violating phase S. 

Because of the weak nature of neutrino interactions 
the use of heavy nuclear targets is unavoidable, and 
this complicates the analysis of the precision physics dis- 
cussed above since model-dependent nuclear corrections 
must be applied to the data. Our present understanding 
of the nuclear corrections is mainly based on charged- 
lepton-nucleus {£A) DIS data. In the early 80s, the Eu- 
ropean Muon Collaboration (EMC) found that the 
nucleon structure functions F2 for iron and deuterium 
differ. This discovery triggered a vast experimental pro- 
gram to investigate the nuclear modifications of the ra- 
tio R[F2^] = Fj"*/(^F|^) for a wide range of nuclear 
targets with atomic number A, see Table ID By now, 
such modifications have been established in a kinematic 



range from relatively small Bjorken x {x ^ 10~^) to large 
X (x ^ 0.8) in the deep-inelastic region with squared 
momentum transfer > 1 GeV^. The behavior of 
the ratio can be divided into four regions: (i) 

R > 1 for X > 0.8 (Fermi motion region), (ii) R < 1 
for 0.25 < a; < 0.8 (EMC region), (in) i? > 1 for 
0-1 ^ 2; ^ 0-25 (antishadowing region), and (iv) R < 1 
for X < 0.1 (shadowing region), with different physics 
mechanisms explaining the nuclear modifications. The 
shadowing suppression at small x occurs due to coherent 
multiple scattering inside the nucleus of a qq pair coming 
from the virtual photon Q with destructive interference 
of the amplitudes 0|. The antishadowing region is theo- 
retically less well understood but might be explained by 
the same mechanism with constructive interference of the 
multiple scattering amplitudes Q or by the application 
of momentum, charge, and/or baryon number sum rules. 
Conversely, the modifications at medium and large x are 
usually explained by nuclear binding and medium effects 
and the Fermi motion of the nucleons 

Instead of trying to address the origin of the nuclear 
effects, the data on nuclear structure functions can be 
analyzed in terms of nuclear PDFs (NPDFs) which are 
modified as compared to the free nucleon PDFs. Relying 
on factorization theorems in the same spirit as in the free 
nucleon case, the advantage of this approach is that the 
universal NPDFs can be used to make predictions for a 
large variety of processes in iA, vA, pA, and AA colli- 
sions. In addition, the nuclear correction factors required 
for the interpretation of the neutrino experiments can be 
calculated in a fiexible way, taking into account the pre- 
cise observable, the nuclear A, and the scale Q^. The 
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Table I. The charged-lepton DIS data sets together with DY 
and with neutrino DIS data sets used in the fit. The table 
details the specific nuclear targets, and the number of data 
points with kinematical cuts (Q^ > 4 GeV'^, W > 3.5 GeV). 
References for the data sets are cited in Refs. 0,01 



ID 


Observable 


A/A' (A) 


ExperimGnt 


# data 


I 




Hc/D 


SLAC-E139, NMC-95 re 


15 


2 


^2 1^2 


Li/D 


NMC-95 


w 


3 


^2 1^2 


Bc/D 


SLAC-E139 


3 


4 


^2 1^2 


C /D 


EMC-88,90, SLAC-E139 
NMC-95, 95 re, FNAL-E665-95 


38 


5 




N /D 


BCDMS-85 


9 




pA i-pA' 


Al/D 


SLAC-E049 E139 


3 


7 


1 F^' 
^2 1^2 


Ca/D 


EMC-90, SLAC-E139 
NMC-95 ic. FNAL-E665-95 


17 


8 


F^ 1 F^' 

^2 1^2 


Fc/D 


BCDMS-85, 87 
SLAC-E049,E139.E140 


24 


9 


F^ IF^' 

-^2 / -^2 


Cu/D 


EMC-88 93 


27 


10 


F^ IF^' 

-^2 / -^2 




SLAC-E139 


2 


11 


^2 I ^2 


Sn/D 


EMC-88 


g 


12 


F^ IF^' 

^2 I ^2 


Xc/D 


FNAL-E665-92 


2 


13 


F^ IF^' 

^2 I ^2 


Au/D 


SLAC-E139 


3 


14 


F^ IF^' 

^2 1^2 


Pb /D 


FNAL-E665-95 


3 


15 


pA ipA' 
-^2 1^2 


Bc/C 


NMC-96 


14 


16 


Fjf'/F^' 
2 12 


Al/C 


NMC-96 


14 


17 


F^ 1 F^' 


Ca/C 


NMC-95, 96 


29 


18 


F^/F^' 


Fc/C 


NMC-95 


14 


19 


F^/F^' 


Pb/C 


NMC-96 


14 


20 


^2 1^2 


C /Li 


NMC-95 


7 


21 




Ca/Li 


NMC-95 


7 


22 


^2 1^2 


He/D 


Hcrrxics 


17 


23 


1 F^' 

^2 1^2 


Kr /D 


Hermes 


12 


24 


F^ 1 F^' 

^2 1^2 


Sn/C 


NMC-96 


111 


25 


F^ 1 F^' 


N/D 


Hcrnics 


19 


32 


F^ 1 F^' 


D 


NMC-97 


201 


26 


VA I P A' 
^ DY 1 DY 


C/D 


FNAL-E772 


9 


27 


pA 1 pA' 
DY 1 DY 


Ca/D 


FNAL-E772 


9 


28 


pA 1 pA' 
" DY / " DY 


Fe/D 


FNAL-E772 


9 


29 


pA j pA 
^ DY 1 ^ DY 


W/D 


FNAL-E772 


9 


30 


pA / pA' 
DY 1 DY 


Fc/Be 


t?M A T T7QC^ 

r iN Aij-rjODD 


28 


31 


pA 1 pA 
^ DY 1 ^ DY 


W/Be 


FNAL-E866 

l^A DIS DY Total: 


28 
708 


33 


da'^^ / dx dy 


Pb 


CHORUS V 


824 


34 


da'^^ /dx dy 


Pb 


CHORUS V 


412 


35 


da^^ 1 dx dy 


Fc 


NuToV V 


1170 


36 


da^^ /dx dy 


Fc 


NuTcV V 


966 


37 


da^^ /dx dy 


Fe 


CCFR di-;i 


44 


38 


da^^ 1 dx dy 


Fe 


NuTeV di-;x 


44 


39 


da^^ 1 dx dy 


Fc 


CCFR di-Ai 


44 


40 


da^^ 1 dx dy 


Fe 


NuToV di-/i 
vA Total: 


42 
3134 



factorization assumption in the nuclear environment is 
therefore a question of considerable theoretical and prac- 
tical importance and global analyses of NPDFs based on 
lA DIS and fixed target Drell-Yan (DY) data confirm its 
validity in the presently explored kinematic range. 

However, in a recent analysis Q of vFe DIS data from 
the NuTeV collaboration we found that the nuclear cor- 



Table II. Summary table of a family of compromise fits. 



w 




(/Pt) 


vA 


(/pt) 


total x^(/pt) 





708 


638 (0.90) 






638 (0.90) 


1/7 


708 


645 (0.91) 


3134 


4710 (1.50) 


5355 (1.39) 


1/2 


708 


680 (0.96) 


3134 


4405 (1.40) 


5085 (1.32) 


1 


708 


736 (1.04) 


3134 


4277 (1.36) 


5014 (1.30) 


oo 






3134 


4192 (1.33) 


4192 (1.33) 



rection factors are surprisingly different from the predic- 
tions based on the l^Fe charged-lepton results with im- 
portant implications for global analyses of proton PDFs. 
This is not completely unexpected since the structure 
functions in charged current (CC) neutrino DIS and neu- 
tral current (NC) electron or muon DIS are distinct ob- 
servables with different parton model expressions. From 
this perspective it is clear that the nuclear correction fac- 
tors will not be exactly the same even for a universal set 
of NPDFs. Note also that some models in the litera- 
ture predict differences between reactions in CC and NC 
DIS [8|. What is, however, unexpected is the degree to 
which the R factors differ between the structure func- 
tions F2^^ and i^^^'^. In particular the lack of evidence 
for shadowing in neutrino scattering down to x ~ 0.02 is 
quite surprising. 

The study in Ref. left open the question, whether 
the neutrino DIS data could be reconciled with the 
charged-lepton DIS data by a better flavor separation 
of the NPDFs. In this letter, we address this question 
in the A-dependent framework of Ref. by performing 
a global x^-analysis of the combined data from lyA DIS, 
M DIS and the DY process listed in Table HI 

When combining neutrino and charged-lepton+DY 
data into a compromise fit, we introduce a weight pa- 
rameter w into the via: 

l^A data i^A data 

The w factor allows us to adjust for the different number 
of points in the separate data sets, and provides a param- 
eter that interpolates between the vA and the £^A+DY 
data. We should stress that the cited in Table llll and 
also in the text is the standard x^; Eq. ((T]) is only used 
internally in the fitting procedure. We construct a set 
of compromise fits with weights w = {0, j, ^, 1, oo} and 
study the dependence of the result on this weight. The 
fit to only neutrino data, denoted w = oo in Table |lll 
is compatible with the results in Q. Similarly, the fit 
to only charged-lepton+DY data, denoted w = 0, agrees 
well with the analysis in Q ■ 

We first compute the nuclear correction factors 
~ Fi^yFi'^ and R[F^^^] ~ F^^^F^'^ in the 
QCD parton model at next-to-leading order employing 
the NPDF fits in Table |lT] for the numerator and free nu- 
cleoli PDFs for the denominator. Q The x-dependence 
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of and R[F^^''] is shown in Fig. □ a) and b), 

respectively, at = 5 GeV^. Similar results hold at 
= 20 GeV^ which we do not present here. We ob- 
serve that the fit to only £A DIS+DY data {w = 0) 
weh describes the SLAG and BCDMS points in Fig. [T] 
a). The same is true for the fit to only DIS data 
(w = oo) which is compatible with the results from the 
NuTeV experiment exemplified in Fig. [1] b) . How- 
ever, comparing the results obtained with the w = and 
the w ~ oo fits one can see that they predict considerably 
different x-shapes. 

The fits with weights w ~ {|, ^, 1} interpolate between 
these two incompatible solutions. As can be seen in Fig.[T] 
a) and b), with increasing weight the description of the 
£Fe data is worsened in favor of a better agreement with 
the I'Fe points. This trend clearly demonstrates that 
the £Fe and the vFe data pull in opposite directions. 
We identify the fits with ■u; = l/2orti; = las the best 
candidates for a possible compromise. To be able to de- 
cisively accept or reject the compromise fits, we apply a 
statistical goodness-of-fit criterion 
probability distribution for the 
N degrees of freedom; 
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Figure 2. Ax^/ACgo as defined in Eq. for the 40 individual 
data sets. Results are shown for the w — ^-fit (a) and tlie fit 
'Ucor5' (b) with w — 1. The solid and dashed lines indicate 
the 90% and 99% confidence limits. The highlighted data 
sets correspond to DIS e^Fe (ID=8), uFe (ID=35), and uFe 
(ID=36). 



llHl3| based on the 



given that the fit has there is only a 10% probability that a fit with > 



(^2)iV/2-lg-xV2 

2^/2r(iV/2) 



(2) 



This allows 



us 

,2 



to define the percentiles via 
/p^" P(X^ Af)dx^ = P% where p = {50,90,99}. Here, 
^50 serves as an estimate of the mean of the x^ distri- 
bution and ^90, for example, gives us the value where 
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Figure 1. Predictions from the compromise fits for the nuclear 
correction factors 7?[F|^"] ~ F|^7F|^ (a) and R[F^^^] ~ 
F^^yp^^ (b) as a function of x for = 5GeV^ The data 
points displayed in (a) are from BCDMS and SLAC experi- 
ments (for references see and those displayed in (b) come 
from the NuTeV experiment [lo| . 



90 

genuinely describes the given set of data. In a global 
PDF fit, the best fit x^ value often deviates from the 
mean value because the data come from different possibly 
incompatible experiments having unidentified, unknown 
errors which are not accounted for in the experimental 
systematic errors. For this reason we rescale the ^90 and 
^99 percentiles relative to the best fit Xo [HI to define 
C90 = Xo(^9o/^5o)and C99 = xli^mUzQ)- This defines 
our criterion: a fit with a given x^ is compatible with 
the best fit with Xo at 90% (99%) confidence if x^ < C90 
(x^ < C'gg)- We apply it to both the total x^ and the 
of the individual data sets. 

For the lA DIS+DY data we use the fit with w — Q a& 
benchmark with Xo = 638 and N — 677 degrees of free- 
dom (for 708 data points and 31 free parameters). The 
upper limits on the x^ at 90% and 99% confidence level 
(G.L.) are then C^g^ = 684 and C^g^ = 722. The bench- 
mark fit for the vA DIS data {w = 00) uses 3134 data 
points with 33 free parameters resulting in iV = 3101 
and one finds C^^ = 4330 and C^^ = 4445. We see 
that none of the compromise fits satisfies both limits at 
the 90% G.L. which is usually used in global analyses of 
PDFs to define the uncertainty bands. At the 99% G.L., 
there are two fits {w = 1/2, w = 1) which are below the 
Cgcf limit. However, only the w = 1/2 fit satisfies the 
corresponding constraint from the charged-lepton bench- 
mark fit. 

We now apply our criterion also to the individual data 
sets with IDs between 1 and 40 in Table H] For the £A 
DIS+DY data (ID=[1,31]) we determine the 31 C90 (C99) 
limits by using the individual xf of the w = fit as Xo i- 
For the vA DIS data (ID=[33,40]) we proceed in a similar 
manner using the individual xf of the u" = 00 fit. The 
results of this detailed analysis are depicted in Fig. [51 
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where we show the quantity 
Ax' 



Xi 



xl. 



AC, 



90 



c, 



90, i 



9 

AO,l 



(*-l,..-,40), 



(3) 



where xi represents the x^-value of the i'th data set. In 
cases where Xi > Cgo^i the fit is not compatible with the 
best fit at the 90% level and Ax^/ACgo > I. The exact 
90% C.L. limit is shown as a constant solid line and the 
dotted line represents the 99% confidence limit. The local 
application of the hypothesis-testing criterion reveals 
that even the compromise fit with weight w = ^ which 
was considered acceptable at the 99% C.L. when looking 
at the nuclear correction factors and at the global change 
in x'^-i cannot be accepted as a compromise solution as 
both the charged-lepton and neutrino DIS data on iron 
exceed the 99% limit. 

In conclusion, the tension between the i^Fe and i/Fe 
data sets leaves us with no possible compromise fit when 
investigating the results in detail, not even when using 
the 99% percentile as the limit as opposed to the more re- 
strictive 90% limit which is usually used to construct the 
error PDFs. This detailed analysis confirms the prelimi- 
nary conclusions of Refs. Q that there is no possible 
compromise fit which adequately describes the neutrino 
DIS data along with the charged-lepton data. 
At face value, this conclusion differs from some results in 
the literature which argue the vA and £^A data are in 
accord Here, we believe an essential element in our 
analysis is the use of the correlated systematic errors of 
the vA data. To highlight this point, we now repeat our 
analysis, but we combine the statistical and all systematic 
errors in quadrature (thereby neglecting the information 
contained in the correlation matrix) for vA data for the 
w = 1 fit with > AGeV (as before); we denote this 
the "Ucor4" fit, and we obtain x' /pt of 1.14 for i^A and 
1.00 for vA. We also use a > bGeV fit (denoted 
"Ucor5") to mimic the cuts of Ref. [l3|; bere we obtain 
x'/pt of 1.14 for l^A and 0.96 for vA. 
If we examine the total x^ values, we find the x' /dof ^ 1, 
and might be tempted to conclude we are able to fit both 
the uA and l^A data simultaneously. However, if we 
look at individual data sets and apply our hypothesis- 
testing criteria, the picture is quite different. Figure [2] 
(b) displays the results for the Ucor5 fit. The higher 
cut of the UcorS fit removes some of the very precise 
NuTeV data at small- a; , thus resulting in an improved x? 
compared to Ucor4. Nevertheless, many of the t'^A data 
sets (ID=3,4,5,6,8) still lie outside the 99% CL percentile. 
Thus, we still conclude that there is no compromise fit 
for the vA and l^A data even if we relax the constraints 
by using uncorrelated errors. 

Consequently, the nuclear correction factor for the neu- 
trino DIS data are indeed incompatible with that of the 
charged-lepton DIS and DY data, and this result depends 
crucially on the use of the precision correlated errors of 



the neutrino data. This result has important implications 
for both nuclear and proton PDFs. If we do not know the 
appropriate nuclear correction to relate different nuclear 
targets, our ability to extract PDFs is limited. For exam- 
ple, the CTEQ6.6 analysis [l^ sidesteps these issues by 
removing most of the i^A data from the fit; however, they 
retain the NuTeV dimuon data since this data is critical 
to constraining the strange quark PDF. This underscores 
the importance of the i/A data for fiavor differentiation. 
Although the NuTeV data provide the tightest con- 
straints due to their statistics, we note that this issue 
cannot be tied to a single data set. For example, we 
find that NuTeV is generally compatible with CCFR and 
CDHSW[ll|. The CHORUS vPh and DPb data have 
larger uncertainties, so they can be compatible with both 
the £^A data and the NuTeV i^Fe data because the 
Ax^/ACgo < 1 for all weights. Compared to the theory 
predictions, NuTeV agrees well in the central x region, 
but exhibits differences both for low x at low Q^, and 
also for very high x {x ^ 0.65). 

We have demonstrated that the vA and £^A data prefer 
different nuclear correction factors, and that there is no 
single "compromise" result that will simultaneously sat- 
isfy both data sets. While we have focused on the phe- 
nomenological aspects for the present study, this result 
has strong implications for the extraction of both nuclear 
and proton PDFs using combined neutrino and charged- 
lepton data sets. Possibilities include unexpectedly large 
higher- twist effects, or even nonuniversal nuclear effects; 
we leave such questions for a future study. 
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